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ABSTRACT 

We study hard states of the black-hole binary XTE J1550-564 during its 2000 
outburst. In order to explain those states at their highest luminosities, L ~ 10% 
of the Eddington luminosity, Le, we propose a specific hot accretion flow model. 
We point out that the highest values of the hard-state L are substantially above 
the L an advection-dominated accretion flow (ADAF) can produce, ~ 0.4« 2 Le, 
which is only ~ (3-4) %Le even for a as high as 0.3. On the other hand, we 
successfully explain the hard states with L ~ (4-10)% using the luminous hot 
accretion flow (LHAF) model. As 10%Le is also roughly the highest luminosity 
an LHAF can produce, such an agreement between the predicted and observed 
highest luminosities provides by itself strong support for this model. Then, we 
study multi- waveband spectral variability during the 2000 outburst. In addition 
to the primary maxima in the optical light curves, secondary maxima were de- 
tected after the transition from the very high state to the hard state. We show 
that the secondary maxima are well modeled by synchrotron emission from a 
jet formed during the state transition. We argue that the absence of the corre- 
sponding secondary peak in the X-ray light curve indicates that the X-ray jet 
emission, regardless of its radiative process, synchrotron or its Comptonization, 
is not important in the hard state compared to the emission from the accretion 
flow. 



Subject headings: accretion, accretion disks — black hole physics — ISM: jets 
and outflows — stars: individual (XTE J1550-564) — X-rays: stars 



1 Shanghai Astronomical Observatory, Chinese Academy of Sciences, 80 Nandan Road, Shanghai 200030, 
China; fyuan@shao.ac.cn 

2 Department of Physics, Purdue University, West Lafayette, IN 47907 

3 Centrum Astronomiczne im. M. Kopernika, Bartycka 18, 00-716 Warszawa, Poland 

4 Department of Astronomy, Peking University, Beijing 100871, China 

5 Joint Institute for Galaxy and Cosmology (JOINGC) of SHAO and USTC 



- 2 - 



1. Introduction 

Black hole X-ray binaries (BHXBs) appear in five main spectral states, namely quiescent, 
low/hard, intermediate, high/soft, and very high ones (e.g., McClintock & Remillard 2006, 
though they prefer a different nomenclature). The origin of the X-ray spectrum of the 
hard state presents a highly interesting problem. The presence of a universal high-energy 
spectral cutoff above ~100-200 keV in the hard state points out to the X-ray emission coming 
from thermal Comptonization in a hot accretion corona with an electron temperature of 
~ 10 9 K ~ 10 2 keV//c (see, e.g., Zdziarski 2000; Zdziarski & Gierlihski 2004 for reviews). We 
note that an alternative jet model (Markoff, Nowak & Wilms 2005), though also relying on 
thermal Comptonization as the dominant component at photon energies > 10 keV, requires 
that the electron temperature is several tens higher, kT ~ 3-4 MeV. Then, the origin of 
the ~100 keV cutoff is not clear; if it is due to the 1st order scattering by those thermal 
electrons, very strong fine tuning is required to achieve the observed cutoff energy. On the 
other hand, higher order scatterings will yield another spectral cutoff at > 3kT ~ 10 MeV, 
which presence remains to be observationally tested. 

There are two main different models of the hot corona. One relies on electron heating by 
magnetic reconnection above the standard thin disk (Liang & Price 1977; Galeev, Rosner, & 
Vaiana 1979). Since the process of magnetic reconnection remains poorly understood, this 
model has not yet been elaborated in detail, although recent numerical simulations have shed 
some light on this problem (Hirose, Krolik & Stone 2006). The other model is a hot accretion 
flow, which is advection dominated (ADAF) in most of its parameter space (Narayan & Yi 
1994, 1995; Abramowicz et al. 1995; Rees et al. 1982). We need to illustrate two points 
here when we use the term ADAF. First, it is found that in most of the parameter space of 
ADAF, outflow is moderately strong, as emphasized by Blandford & Begelman (1999, the 
ADIOS solution) and other authors. Second, as we will show later, another hot accretion 
solution may be responsible for the X-ray emission of the relatively luminous hard states. In 
comparison with the magnetic reconnection, ADAF has clear dynamics, and its astrophysical 
applications can be worked out in detail. The specific suggestion of applying the ADAF to 
the hard state is due to Narayan (1996) and Narayan, McClintock & Yi (1996) (see also 
Shapiro, Lightman & Eardley 1976, who, however, did not recognize the importance of 
advection). First detailed calculations applied to data were then done by Esin, McClintock 
& Narayan (1997, hereafter E97). 

In the present paper we concentrate on the hot accretion flow. We note that the radiative 
efficiency of an ADAF increases with the accretion rate, M. Thus an ADAF can explain not 
only very dim sources like the Galactic center (e.g., Yuan, Quataert & Narayan 2003), but 
also some relatively luminous sources such as hard states of BHXBs. The highest rate, called 



- 3- 



the critical accretion rate, of an ADAF is determined by the balance between the viscous 
heating and cooling, g vis « q ie , which occurs at Madaf ~ 10a 2 M E , where a is the viscosity 
parameter, Me = Le/c 2 , q is the rate of energy change per unit volume, and Le is the 
Eddington luminosity. The corresponding maximum luminosity is Ladaf ~ 

0Aa 2 L E (E97)H 

The actual value of a remains uncertain, but recent three-dimensional magnetohydrodynamic 
(MHD) numerical simulations suggest a ~ 10 _2 -10 _1 in most of the disk midplane (Hawley 
& Krolik 2001). Even for a as high as 0.3 (a common value used in ADAF studies), the 
maximum ADAF luminosity is still only ~ (3-4) %L E . It has been thought that above 
the Ladaf there is no hot accretion solution, with the thin disk remaining the only viable 
solution. 

However, observations of the hard state in BHXBs often show higher luminosities during 
the rising part of an outburst. For example, the hard state of XTE J1550-564 reached 20%Lg 
(Sobczak et al. 2000) and ~ 10%Le (Rodriguez, Corbel & Tomsick 2003; this work) during 
its outbursts in 1998 and 2000, respectively. The highest hard-state luminosity of GX 339-4 
was ~ 25-30%Le (Zdziarski et al. 2004; Done & Gierlihski 2003). These observational 
results have long presented a puzzle for the hot accretion flow models. 

Of importance in this context was the finding of a new hot accretion flow solution, the 
so-called luminous hot accretion flow (LHAF; Yuan 2001, hereafter Y01). Different from 
an ADAF, this solution corresponds to accretion rates above the critical rate of ADAF, 
-M\daf j up to another critical accretion rate which is about 3-5 times higher than Madaf 
for parameters adopted in Y01. In addition to accretion rates, the radiative efficiency of an 
LHAF is also higher than that of a typical ADAF. Therefore, an LHAF can produce much 
higher luminosity than an ADAF. Yuan & Zdziarski (2004) studied some luminous hard 
states of BHXBs and AGNs and found that their relatively low electron temperatures and 
high luminosities do require an LHAF to be present. In the first part of this work, we apply 
this approach to hard-state spectra of XTE J1550-564 during its 2000 outburst. 

In addition to studying X-ray spectra alone, more and more efforts are devoted to 
multi- wavelength simultaneous observations and modeling of BHXBs, as illustrated, e.g., by 
studies of XTE J1118+480. Its nearly simultaneous observations were conducted from radio 
to X-ray wavelengths (e.g., Hynes et al. 2000; McClintock et al. 2001; Frontera et al. 2001, 
2003; Chaty et al. 2003). These have allowed important theoretical conclusions to be drawn. 
Esin et al. (2001) showed that the EUV data rule out models with a thin disk extending to 
the last stable orbit. Using a coupled jet-ADAF model, Yuan, Cui & Narayan (2005) and 



lr Thus the highest radiative efficiency of an ADAF is ~ 0.04, not much lower than that of a standard thin 
disk around a Schwarzschild black hole, ~ 0.057. 
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Malzac, Merloni & Fabian (2004) successfully explained the radio to X-ray spectrum and 
most of the timing properties of the source. 

Here, we apply the same approach to XTE J1550-564. We use simultaneous multi- 
wavelength observations at radio, infrared, optical and X-rays during the 2000 outburst 
(Corbel et al. 2001; Jain et al. 2001, hereafter J01; Tomsick, Corbel & Kaaret 2001; Rodriguez 
et al. 2003). Compared to XTE J1118+480, XTE J1550-564 is much more luminous. During 
the outburst, transitions between the hard and very high states occurred, unlike the case 
of XTE J1118+480, which was always in the hard state. Therefore some new interesting 
phenomena were observed. One issue we discuss below is that of secondary maxima in the 
infrared/optical light curves. J01 reported infrared (H), optical (/, V) and X-ray light curves 
for the 2000 outburst. After reaching the peak of the outburst, the IR, optical, and X-ray 
fluxes declined, and a transition from the very high state to the hard state occurred around 
MJD 51,680. Afterwards, the X-ray flux continued to decline whereas the IR/optical light 
curves began to rise and reached secondary maxima at MJD 51,699-51,705, after that they 
declined again. A contemporaneous spectrum in the radio band was obtained by Corbel et 
al. (2001) at MJD 51,697. Here, we model the simultaneous radio, IR, optical, and X-ray 
data and explain the secondary IR/optical maxima. 

In $21 we model the X-ray emission of the hard states, paying special attention to the 
dynamics of the hot accretion flow. In §3j we study the simultaneous radio, IR/optical and 
X-ray data during the decline phase and explain the secondary IR/optical maxima. We 
conclude in §5]with a summary and discussion. 

2. Modeling the X-ray spectra of the hard states of XTE J1550-564 

2.1. Observational data 

During its 2000 outburst, XTE J 1550-564 experienced hard and very high states. We 
have selected four occurences of the hard states, at MJD 51,646, 51,658, 51,687, and 51,696, 
based on the criteria that they span a wide range of X-ray luminosity and that each includes 
simultaneous optical data. Their spectra are shown in Figure [TJ by different colors. For 
comparison, we also show an example of the very high state spectrum, of MJD 51,662. Their 
X-ray RXTE spectra were earlier studied by Rodriguez et al. (2003) and Tomsick et al. 
(2001); here we have reextracted them ourselves. They can be well fitted by a power-law form 
with an exponential cut-off. The X-ray spectra in blue, red, green, and magenta in Figured] 
(from top to bottom) are fitted with the photon index and the e-folding energy, (T,E{), of 
(1.70 ± 0.01, 115 ± 5.6), (1.46 ± 0.01, 1371?), (1-53 ± 0.01, 200ljj{j) and (1.53 ± 0.01, 460l^o)> 
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respectively. We see a strong anticorrelation between the X-ray flux and Ef, with its value 
monotonously increasing with the decreasing flux. 

The IR/optical data in Figured] are from J01. Since we do not know the exact value of 
the extinction, we show two sets of fluxes for each spectrum, corresponding to the two limits 
of the extinction of Ay = 4.75 (Orosz et al. 2001) and Ay = 2.2 (Sanchez- Fernandez et al. 
1999). The details of the data reduction can be found in Xue, Wu, & Cui (2007). 



2.2. Accretion flow models 

We now fit the spectra using our accretion flow model. The components of it are nearly 
the same as in E97, i.e., the accretion flow consists of an inner hot accretion flow (ADAF or 
LHAF) within a transition radius, r tr , and an outer thin disk. We assume that only a fraction 
of the accretion rate at r tr actually accretes onto the black hole and the rest is ejected. This 
assumption is required by both numerical simulations (Stone, Pringle, & Begelman 1999; 
Hawley & Balbus 2002; Igumenshchev et al. 2003) and analytical work (Narayan & Yi 1994; 
Blandford & Begelman 1999; Narayan et al. 2000; Quataert & Gruzinov 2000). Details of 
this process depend on the accretion rate. Here, in order to formally take into account the 
role of outflows as well as convection in modifying the density profile of the accretion flow, 
we assume following Yuan et al. (2005), that M(r) = M for r > r tr , and 

d In M(r) 

dlnr = s ( r )> r < r tn (!) 

where 

s(r) = s max[/(r),0], (2) 

and So is independent of r but it can be different for different values of M . The above 
formula is based on the consideration that the outflow (and convection) is ultimately due to 
the accreting gas acquiring a positive Bernoulli parameter, as emphasized by Narayan & Yi 
(1994), and which value is smaller for higher accretion rates due to the stronger radiative 
energy loss. Then, /(r) is the advection factor of the accretion flow, defined as, 

f(r) = — = — • (3) 

yVis yVis 

A negative value of / in eq. ([3]) implies that advection plays a heating rather than a cooling 
role. In this case, the hot accretion flow is in the LHAF regime (Y01), as we will point out 
in §2.3. When M is very low, f(r) = 1 and s(r) = s - m this case, eq. (JTJ) gives us the usual 
form, M = Mo(r/r tr ) s ° (e.g., Blandford & Begelman 1999). 
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We calculate global solutions of the hot accretion flow using the method of Y01. A major 
difference of that method with respect to that of E97 is that we solve the radiation hydro- 
dynamics equations self-consistently, and thus obtain f(r) at each radius. In contrast, E97 
used the approximation of f(r) having a constant average value at all radii. The radiation 
processes we consider include bremsstrahlung, synchrotron emission, and Comptonization of 
both synchrotron photons from the hot accretion flow and soft photons from the thin disk. 
We use formulae in Narayan & Yi (1995) and Coppi & Blandford (1990) to calculate the 
synchrotron and Comptonization spectra, respectively, as a function of the radius. We treat 
self-absorption as in Manmoto, Mineshige & Kusunose (1997). The emission from the outer 
cool disk is modeled as a multicolor blackbody spectrum. The temperature as a function 
of r is determined by the viscous dissipation and the irradiation of the disk by the inner 
hot flow. E97 show that when r tr > 30rs, the outer thin disk has little effect to the X-ray 
emission from the inner hot accretion flow because the seed photons for the Comptonization 
mainly come from the synchrotron and bremsstrahlung emissions in the hot accretion flow. 
Here r$ = 2GM/c 2 is the Schwarzschild radius. 

Now we discuss the values of the model parameters. As we stated in §TJ the exact value 
of a is uncertain, and we adopt a = 0.3, which appears the highest possible value. This 
choice allows us to show the highest values of L possible with the ADAF model, and it shows 
that even then that model cannot explain the luminous hard states of BHXBs. We set the 
magnetic parameter, (3 defined as the ratio of the gas pressure to the sum of the gas and 
magnetic pressures, at /3 — 0.9 (from MHD numerical simulations; e.g., Hawley & Krolik 
2001), and the fraction of the viscous dissipation directly heating electrons is set at 5 = 0.5 
(from the modeling to the supermassive black holein our Galactic center: Yuan et al. 2003). 
The parameters a, (3 and S are not free hereafter, though we realize that their values bear 
large uncertainties. The outer radius of the cold disk is set equal to the tidal radius, which 
is ~0.87 of the Roche lobe radius (Papaloizou & Pringle 1977). The Roche lobe radius is 
0.55-0.59 of the separation (Eggleton 1983) for the estimated mass ratio of the black hole 
to the companion of 6.7-11, while the separation is ~ 9 x 10 11 cm (Orosz et al. 2002). This 
implies the outer radius ~5x 10 11 cm, equivalent to ~ 1.5 x 1O 5 (1OM /M)rg. We set the 
black-hole mass to M = 1O.5M and the distance to D = 5.3 kpc (Orosz et al. 2002)0. 

Another important parameter is the transition radius, r tr . As it has been shown by 
Yuan et al. (2005), its value is best constrained by the thin disk emission close to r tr , which 
is in the EUV/soft X-ray band. However, we have no data for that band, and hereafter 



2 Jonker & Nelemans (2004) argued that the distance to XTE J1550-564 might have been underestimated 
by as much as a factor of three due to an overestimation of the interstellar extinction. However, if it were 
true, the highest 1-20 keV luminosity of the 1998 outburst would be ~ IOLe, which seems highly unlikely. 
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assume r tr = lOOrs- As stated above, most of the seed photons for Comptonization at this 
value of r tr are synchrotron photons from the hot accretion flow. Our satisfactory fits to 
the data (the magenta, green and red colors in Figure [T]) indicate that our solution is viable 
though it may be not unique. 

On the other hand, it is possible that when the hard state approaches the very high state, 
the soft photons from the thin disk or clumps in the hot accretion flow become dominant 
(see also Wardzihski & Zdziarski 2000). As it has been shown by Yuan et al. (2007), the 
X-ray slope of (at least some) black-hole binaries is correlated with the luminosity, L, and 
the sign of the correlation changes when the luminosity of the hard state increases above a 
certain value. Below this luminosity, the X-ray spectrum hardens with the increasing L, and 
above it, it softens. The boundary in our case approximately corresponds to the red-color 
spectrum in Figure dj though in our case with a limited range of L the X-ray spectra below 
it keep an approximately constant slope. The change of the sign of the correlation may be 
due to the change of the dominant source of the seed photons. 



2.3. Modeling results 

Figure [T] shows our X-ray spectral fits to the lower three hard states. (We will discuss 
the brightest hard state separately.) From bottom to top, they are shown by the dot-dashed, 
dashed and solid curves and have (M /Me,so) equal to (1.0,0.55), (1.1,0.3), and (1.3,0.3), 
respectively. The advection factor, /, and the electron temperature, T e , of these models are 
shown as a function of r in Figures [2] and [3j respectively. 

We see in Figure [TJ that the three models fit the X-ray spectra relatively well, including 
both the spectral slope and the cut-off energy. E97 did not take into account outflows and 
predicted that the X-ray spectra from hot accretion flows should harden with the increasing 
accretion rates. The X-ray spectrum at MJD 51,646 (in red) is harder than that at MJD 
51,687 (in green), consistent with that prediction. However, the slope of the spectrum at 
MJD 51,687 is identical to that at MJD 51,696 (magenta), different from the prediction of 
E97. Our calculations show that the inclusion of the outflow in our model can mostly resolve 
this discrepancy. From eqs. ([IH2]) and Figure [21 we see that there is no outflow for the model 
shown by the (high) solid curve since its advection factor is / < 0. For the (middle) dashed 
line model, / is very small: / ~ 0, so outflow is also weak. Since the accretion rate of the 
solid curve is larger than that of the dashed curve, the former spectrum is harder. On the 
other hand, for the (low) dot-dashed curve, so = 0.55 and / ~ 1, the outflow is strong. 

Then, the presence of an outflow will make the emitted spectrum harder compared to 
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that with no outflow at the same Mq. For comparison, we also show in Figure[T]the spectrum 
produced by a model without outflow (at M = 0.5M E to reproduce the X-ray flux level), see 
the double-dot-dashed curve. We clearly see that its spectrum is softer than that of the dot- 
dashed curve (including the outflow). The physical reason for this is as follows. The emitted 
X-ray spectrum is the sum of the local Comptonization spectra at each radius, and the slope 
of that spectrum is determined by the Compton y parameter, with a larger y corresponding 
to a harder spectrum (e.g., Beloborodov 1999). At small radii, y is relatively small so the 
corresponding spectrum is soft. When a strong outflow exists, the accretion rate decreases 
with decreasing radii, thus the contribution of the innermost region (where the dissipation 
peaks) becomes smaller. Therefore, the total spectrum will become harder. We note that 
this effect is not universal; it holds only for a certain range of the accretion rate. 

We note that the spectra produced by all the three models are somewhat too soft in 
comparison with the X-ray data. Increasing slightly Mq can make the spectrum harder, thus 
improving the fits. However, the corresponding luminosities will be higher than observed 
values. Given the uncertainty in the distance of XTE J1550-564, it is possible that the 
distance is slightly larger than D = 5.3 kpc. Another solution would be to use a larger 
(3, which will reduce the magnetic field and thus the synchrotron emission. As the latter 
provides seed photons for Comptonization, this will lead to hardening of the spectra. 

Our models also fit satisfactorily the e-folding energies, Ef. As noted in §2.1[ , Ef is 
anticorrelated with the X-ray luminosity. Such an anticorrelation appears common in the 
hard states. It takes place, e.g., in GRO J0422+32 (Grove et al. 1998) and GX 339-4 
(Wardzihski et al. 2002). In thermal Comptonization models, Ef is close to the electron 
temperature T e . Such an anticorrelation between T e and L (or M) is expected in the hot 
accretion flow models, both ADAF and LHAF, as shown here by Figure [3] (see also Fig. 3b of 
E97). It shows the T e profiles of the three models shown in Figured], together with the values 
of Ef of the observed X-ray spectra. Note the value of Ef for the lowest considered state is 
not precisely determined. Since the emitted X-ray spectra are integrated over the radii, the 
Ef here only measures an average T e . Still, the shown agreement is excellent, indicating that 
the hot accretion flow model is able to quantitatively predict the correct Ef. We emphasize 
that this result is quite robust, independent of the model parameters. It provides strong 
support for this model of the hard state. 

Let us now discuss the dynamics of the accretion flow. We first need to illustrate the 
physics of the LHAF. The energy equation of ions in hot accretion flows reads: g a( j v = 
<?int — <?com P = <?vis — <?ic- The subscripts denote the energy advection, the gradient of the 
internal energy of ions, compression work, viscous heating, and Coulomb collision cooling, 
respectively. The critical accretion rate of an ADAF Madaf is determined by g v i S = g; e 
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(Narayan, Mahadevan & Quataert 1998). Below Madaf, <?vis > Qie, therefore we have qi nt = 
?com P + Qvis — lie > 0. The positive sign of q mt means that the flow can remain hot if it starts 
out hot (because q int oc vdT/dr, and v < 0). This is why a hot ADAF solution can exist. 
Above Madaf up to another (higher) critical accretion rate determined by q vis + g C om P = Qie, 
the value of q- m t(— q C om P + <?vis — <?ie) is still positive. Thus a new hot accretion solution, i.e., 
LHAF, exists between Madaf and this new rate (Y01). For an ADAF, g a d v > 0, so advection 
plays a cooling role. While in an LHAF, g adv < 0, so advection plays a heating role. In this 
case, the radiative cooling is balanced by the sum of advective and viscous heating. 

In FigureEl we see that the advection factor for the (low) dot-dashed curve is > 0, thus it 
represents an ADAF. The bolometric luminosity of this model is L ~ 0.8%Le- The advection 
factor for the (middle) dashed curve is negative in some region, and thus it is marginally an 
LHAF, with L ~ 3%Le- This result is roughly consistent with those of E97, who found the 
highest accretion rate of ADAFs as Madaf ~ 13« 2 Me ~ Me and the corresponding highest 
L is ~ 0Aa 2 L E « (3-4)%L E for a = 0.3. Thus, for spectra such as that of MJD 51,687, 
whose L is ~ 3%Le, an ADAF marginally does not work. For the spectrum at MJD 51,646, 
the luminosity is still higher, ~ 6%Le- In this case, the flow is an LHAF in the entire range, 
see the solid curve in Figure [21 

However, even an LHAF cannot explain the spectrum at MJD 51,658 (the blue sym- 
bols). This is because the required L is above the maximum accretion rate of LHAF, Mlhaf, 
roughly determined by g vis + g com p = fe- Because g com p does not depend on a, the value 
of Mlhaf depends only weakly on a, unlike Madaf- In addition to a, other model pa- 
rameters such as (3, 5 and the outer boundary conditions also affect its value. Y01 found 
that Mlhaf/Me — 5 and 3 for a = 0.3 and 0.1, respectively. For the parameters adopted 
here, Mlhaf is slightly larger than 1.3MeB 

which corresponds to L < 8%L E . This is not 
high enough to explain the X-ray spectrum at MJD 51,658. When M > Mlhaf, the inner- 
most region of the LHAF will collapse and form an cold annulus (Y01). Calculations of the 
emitted spectrum in this case are complicated; here, we hypothesize that this effect may be 
responsible for the very high state. 

Phenomenologically, the X-ray spectrum at MJD 51,658 appears to be intermediate 
between the hard state spectra, and the very high state spectrum, shown by the black 
points. The approximate agreement between the observed highest luminosity of the proper 
hard state and the highest LHAF luminosity strongly supports the LHAF model for the hard 



3 As a comparison, the critical rate of the ADAF for the parameters adopted in the present paper is 
Madaf — 1-1-Me. The Mlhaf /Madaf ratio is smaller than that obtained in Y01. This is due to the 
different parameters assumed, a, (3 and S as well due to the inclusion of outflows in the ADAF but not in 
the LHAF (where they are suppressed). 
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state of BHXBs. 

Figure [T] shows that the IR/optical fluxes are under-predicted by the hot flow models. 
This implies that some other components such as a jet and the thin disk also contribute. We 
consider the role of the jet in §31 Here, we show the emission from an irradiated truncated 
thin disk with M = 1.3Me and r tr = lOOrs. We see that this improves the agreement of the 
model with the data in the optical band, but most of the IR data are still not explained. 

3. The role of jet 

We note here two interesting results in the IR/optical and X-ray light curves for the 
2000 outburst reported by J01. The first one is that the outburst in the IR/optical started 
some days before that in X-rays (seen by the RXTE/ ASM). Here, we focus on another 
interesting result of J01. After the primary peak of the outburst, the HIV and X-ray fluxed 
declined, and the transition from the very high state back to the hard state occurred around 
MJD 51,680. After the transition, the X-ray flux continued to decline while the HIV fluxes 
began to rise and reached secondary maxima at MJD 51,699-51,705. J01 also found that 
the secondary maximum is most prominent in the H band. 

The IR/optical and X-ray spectra at MJD 51,696 are shown in Figure HI It also shows 
the contemporaneous radio measurements at MJD 51,697 (Corbel et al. 2001) with the 
Australia Telescope Compact Array (ATCA) at four frequencies from 1384 to 8640 MHz. 
The detected spectrum is slightly inverted with the spectral index of «r = 0.37 ±0.10 (where 
the energy flux is F v cx z/ QR ). Since the radio regime is much less variable than the X-rays, 
we assume that the radio spectra at MJD 51,696 and 51,697 were the same. 

Based on Figure HI we can rule out the standard thin disk as the origin of the secondary 
HIV peaks. The flux produced by the thin disk in the H band depends mainly on the values 
of the outer radius of the thin disk, r out , and M. The dot-dashed curve in Figure H] shows the 
emitted spectrum by an irradiated thin disk with Mq = I.OMe and r out = 1.5 x 10 5 rs- We 
can see that the emission from the disk (in the Rayleigh- Jeans regime) is substantially below 
the observed H flux. We note that only for M ~ 40Me, which appears highly unlikely, the 
observed and predicted H fluxes would become equal. 

Jet formation is usually observed when the source enters the hard state (Fender 2001, 
2006). Prompted by this, J01 proposed that the secondary peaks are due to the emission 
from the jet formed at the time of the very high-hard state transition around MJD 51,680. 
Our detailed calculations support this conjecture. We model the MJD 51,696 spectrum 
(Figure HJ) taking into account the contributions from the hot flow, the truncated disk, and 
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the jet. We use the internal shock scenario (widely adopted in the study of GRB afterglows) 
to calculate the jet emission, see Yuan et al. (2005) for details. Briefly, internal shocks 
within the jet occur due to collisions of shells with different velocities. The shocks accelerate 
a fraction of the electrons into a nonthermal, power-law, energy distribution. The steady- 
state electron energy distribution is then self-consistently determined taking into account 
the radiative cooling. We parametrize the model by two parameters, e e and 6b, giving the 
fraction of the shock energy going into the accelerated electrons and the (amplified) magnetic 
field, respectively. We then calculate the nonthermal synchrotron emission. 

The modeled jet and total spectra are shown by the dashed and solid curves respectively 
in Figure HI We can see that this model fits the data from the radio to the X-rays very well. 
The jet parameters are: e e = 0.06, e B = 0.02, its opening angle of 9 = 0.1, the bulk 
Lorentz factor, Fj et = 1.2, and the spectral index at which the electrons are accelerated of 
p = 2.23. The values of the parameters are the same as those used by Yuan et al. (2005) 
to model the outburst hard state of XTE Jll 18+480. The mass flow through the jet is 
Mj et = 3.3 x 10 _3 A / /e. This is ~ 0.6% of the accretion rate at in the inner flow, equal to 

0. 54Me at 5rs- This is very similar to the corresponding value of 0.5% in XTE J1118+480. 

Based on the above fitting, we now explain the secondary HIV peaks as follows. The 
jet began to form at the very-high to hard state transition, at MJD 51,680. Initially, M- ]et 
was very small, therefore the contribution of the jet emission to the H flux was negligible. 
But then the accretion M continued to decrease and the accretion luminosity declined. On 
the other hand, the jet grew, with Mj et increasing, so the relative jet contribution to the 
HIV bands continued to increase. The increase of Mj et lasted for ~20 d and reached its 
maximum at MJD 51,699-51,705, which corresponds to the secondary peaks in the HIV 
light curves. After that time, Mj et decreased together with the M. Our fitting results at 
MJD 51,696, shown in Figure HJ are for the time very close to that of the secondary HIV 
peaks. 

Fig. 1 in J01 shows that the H band fluxes at MJD 51,680 and 51,696 differ by one 
magnitude. We note that the modeled ADAF and the thin disk fluxes could have changed 
only very little between those two days (see Figured]). On the other hand, from Figure HI 
we see that the H band flux increases due to the appearance of the jet emission at MJD 
51,696 by about one magnitude, completely consistent with the observational results in J01. 
We also see in Figure H] that while the vF v spectrum from the jet is almost flat in the HIV 
band, the spectrum of the accretion flow (the thin disk and ADAF) rises with the increasing 
frequency. This explains why the secondary peak is most prominent at the H band (J01), 

1. e., at the lowest energy. 

Another important observational result is the absence of the corresponding secondary 
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peak in the X-ray light curve. This requires that whatever emission the jet produces in 
X-rays, it has to be negligible compared to that from the hot accretion flow. Our theoretical 
modeling is consistent with that result. As shown in Figure HI the contribution of the 
jet emission in the X-ray band is much weaker than that of the accretion flow. We thus 
conclude that the absence of the secondary maximum in the X-ray light curve provides a 
strong evidence that the jet emission (synchrotron plus Compton) is not important in the 
X-ray band in the hard state. 



4. Summary and Discussion 

We have investigated the hard states of XTE J1550-564 during its 2000 outburst. We 
have modeled the X-ray spectra with accretion flow models O. Our model consists of an 
inner hot flow within a transition radius and an outer thin disk, as in E97. Different from 
E97, we have taken into account two new developments of the hot accretion flow models. 
One is the inclusion of outflow, the other is the extension of the ADAF model to higher 
accretion rates using the LHAF model (Y01). Our model fits well the lower three of the four 
characteristic hard-state spectra, see Figure HJ The two extensions of the model are necessary 
for the obtained good fits. First, the inclusion of the outflow allows us to fit the X-ray spectral 
slopes. Without outflow, the emitted spectrum hardens with the increasing accretion rate 
(E97), which disagrees with the observations, where the slopes of the two lowest data sets 
are almost identical in spite of their different luminosities. To resolve this discrepancy, the 
outflow becomes weaker with the increasing accretion rate. This is physically motivated as 
the radiative energy loss is smaller at lower accretion rates and thus the Bernoulli parameter 
is larger. 

The extension of hot accretion flow models from ADAFs to LHAFs enables us to model 
the luminous hard states with luminosities up to ~ 8%Le (while the ADAF is limited to 
< 3%L E ). The LHAF fits the spectral slope and the cut-off energy of those states, as shown 
in Figure dj confirming the prediction of Yuan & Zdziarski (2004). 

The value of the highest possible luminosity of the LHAF model, Llhaf, depends on a 
and other model parameters, and the outer boundary condition (see Y01). For our assumed 
parameters, we find Llhaf — 8%Le, somewhat lower than L of the most luminous hard 
state considered here shown by the blue symbols in Figure [TJ However, we can see that this 
state appears intermediate between the hard and the very high states, and perhaps some 
prcesses responsible for the very high state are already for importance here. Thus, the overall 
agreement between the maximum observed L and the prediction of the LHAF model provides 
strong support for it as a model of the X-ray emission of the hard states of BHXBs. 
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The highest L of the hard state varies among various sources as well as different outbursts 
of a given source. As discussed in SJH the maximum hard-state L of XTE J 1550-564 and 
GX 339-4 varies in the ~10-30% range. On the other hand, the corresponding value in Cyg 
X-l is only ~ 0.01-0.02L E (d/2 kpc) 2 (10M©/M) (e.g., Zdziarski et al. 2002). We discuss now 
possible reasons for this variability. 

First, L L haf is a function of the model parameters, e.g., a, and parameter variations 
among different sources will contribute to its dispersion. However, numerical simulations 
of accretion flow indicate that a < 0.3 (e.g., Hawley & Krolik 2001). Thus, this effect can 
explain a low value of the maximum L, but not values > 10%Le- Still, we can obtain 
~ 15%Le by adjusting the outer boundary condition (determined by the physics of the 
transition region between the hot and cool accretion flows, which is complex and currently 
unclear) . 

An efficient way of increasing the maximum hard-state L can be due to the angular 
momentum of the black hole being > 0. In the present calculations, we have assumed a 
Schwarzschild black hole. If the black hole is rapidly rotating, substantially more gravita- 
tional energy will be released because of the deeper potential well (e.g., Gammie & Popham 
1998; Popham & Gammie 1998). Shafee et al. (2006) indeed have found that a number of 
BHXBs appear to have the spin parameter of ~0.6-0.8. 

Yet another effect that could increase I/lhaf is a magnetic coupling between the rotating 
black hole and the accretion flow. If the black hole rotates faster than the disk, energy 
and angular momentum will be extracted from the black hole and transferred to the disk 
by the magnetic field (Li 2000, 2002; Wang et al. 2003). However, we are not able to 
estimate quantitatively the importance of this effect here because the existing studies have 
only considered the standard thin disk rather than a hot flow3. 

As we point out above, the most luminous hard states border the very high state. 
The nature of the latter remains an open question. Numerical calculations by Y01 (see 
also Begelman, Sikora & Rees 1987) indicate that above Mlhaf the hot accretion flow will 
collapse within a certain radius and form a cold annulus (see Liu, Meyer & Meyer- Hofmeister 
2006 for a similar result in the context of their evaporation model). The inner cold annulus 
will emit black body radiation which may explain the thermal component of the very high 
state. At the collapse region, strong magnetic reconnection events should occur due to the 



4 Recently attention has been paid to the effect of the presence of magnetic stresses near and inside the 
last stable circular orbit on the radiative efficiency of the standard thin disk (e.g., Agol & Krolik 2000; 
Beckwith, Hawley & Krolik 2006). This effect does not apply to the hot accretion flow since the no-torque 
boundary condition is not adopted in the global solution of hot accretion flow. 
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overflow of the magnetic field (Begelman et al. 1987; Y01). This is likely to accelerate some 
electrons into a power-law distribution. Compton scattering of the soft photons from the 
inner annulus by both the thermal and power-law electrons may then be responsible for the 
X-ray emission of the very high state. We note that Gierlihski & Done (2003) have found 
that the shape of the soft 7-ray spectrum of XTE J1550-564 in the very high state indeed 
requires the presence of both thermal and nonthermal (power-law like) electrons. 

A possible problem with the above scenario is that the hot accretion flow is strongly 
thermally unstable at the collapse region (Yuan 2003). As a result of the instability, the 
accretion flow may form a two-phase flow structure, with cold clumps surrounded by the 
hot phase (Yuan 2003; Guilbert & Rees 1988; Celotti, Fabian & Rees 1992; Kuncic, Celloti 
& Rees 1997; Krolik 1998; Beloborodov 1999; Merloni et al. 2006). The radiation from 
the cold clumps can then explain the thermal component of the very high state, whereas 
Comptonization of the soft photons by the high-energy electrons in the hot phase may be 
responsible for the power law tail. Thus, this represents an alternative scenario for the very 
high state. (However, a hybrid electron distribution in the hot phase is required as stated 
above, see Gierlihski & Done 2003.) We speculate that even some luminous hard-state 
flows with M < Mlhaf may also contain such a two-phase configuration. The formation of 
the cold clumps in this case could be due to the strong instability in the transition region 
between the hot accretion flow and the standard thin disk (Kato 1999, 2000). This scenario 
is promising in explaining the observed power spectra of the rapid aperiodic variability of X- 
ray emission of BHXBs (Bottcher & Liang 1999). It also may explain the claimed relativistic 
iron K emission line profiles in BHXBs (e.g., Miller et al. 2006). We plan to investigate these 
issues in future work. 

We have also investigated the nature of the secondary peaks in the HIV light curves of 
XTE J1550-564 (J01). They began to emerge after the transition from the very high state 
back to the hard state, and reached the maximum ~20 d later. However, there have been 
no corresponding peak in the X-ray light curve. Based on a detailed modeling to the multi- 
waveband spectrum from the radio to X-rays, we conclude that the secondary peaks are due 
to the emission from the jet formed during the state transition. The mass loss rate in the jet 
kept increasing for ~20 d after its birth, and then decreased. The secondary peaks were then 
due to the corresponding maximum of the relative contribution of the jet {HIV) emission 
with respect to the accretion emission. The absence of the corresponding peak in the X-rays 
strongly indicates that the jet emission (either synchrotron or Compton) is negligible in the 
X-ray band in the hard state. This differs from the scenarios with the dominance of the jet 
in X-rays (e.g., Markoff et al. 2005). In principle, it is possible that a drop in the X-ray flux 
from the disk is exactly compensated by a rise in the X-ray flux from the jet, with the total 
flux changing smoothly. As this scenario requires fine tuning, we consider it to be highly 
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unlikely. 

We also point out that Xue et al. (2007) present another argument against the jet 
dominance in X-rays in the hard state. By modeling the spectrum of the jet blob, they 
obtain the parameters of the electron energy distribution. Assuming the electron energy 
distribution is the same for all blobs in the jet, they model the overall radio spectrum of the 
source and extrapolate their fit to higher energies. They find that the synchrotron emission 
from the jet contribute negligibly to the observed X-ray emission. 

We thank John Tomsick and the referee, Chris Done, for valuable comments. This work 
was supported in part by the National Natural Science Foundation of China (grant 10543003) 
and One-Hundred- Talent Program of China (FY), and the Polish grants 1P03D01827, 1P03D01128 
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Fig. 1. — Spectral fitting results for the hard states of XTE J1550-564 during its 2000 
outburst at different dates. The two sets of optical data for each date correspond to the 
two limiting values of the extinction, Ay = 4.75, 2.2. The parameters of the hot accretion 
flow models are: M = I.OMe and So = 0.55 (dot-dashed curve); M = I.IMe and so = 0.3 
(dashed curve); and M = 1.3Me and sq = 0.3 (solid curve). The dot-dashed curve shows an 
ADAF model while the other two models are LHAFs (see Figure [2]) • The double-dot-dashed 
curve shows a model without an outflow to illustrate its hardening effect on the X-ray slope 
(see §2.31) . The dotted curve shows the emission from a truncated thin disk at M — 1.3 Me 
and r tr = lOOrs- We see that the inclusion of this component imroves the fit in the optical 
range, but it is not sufficient to explain the IR data. 
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Fig. 2. — The profile of the advection factor (defined in eq. EJ) for the three models shown 
in Figure [TJ The model denoted by the dot-dashed curve is an ADAF since its advection 
factor is always positive, while the other two models are LHAFs since their advection factor 
are negative in some regions. The irregular behavior near lOOrs is due to an effect of the 
boundary condition. 
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Fig. 3. — The profile of the electron temperature for the three models shown in Figure [TJ 
With the increase of the accretion rates, the electron temperature decreases, consistent with 
the observed anticorrelation between the e- folding energy, Ef, and the luminosity seen in 
Figure [TJ The three circles with different colors show the Ef for the corresponding three 
X-ray spectra shown in Figure [TJ Their abscissa is arbitrary. The magenta circle (open) 
denotes the e-folding energy poorly constrained by the data. Since Ef is proportional to 
the electron temperature, these results show that the hot accretion flow models (ADAF and 
LHAF) can quantitatively predict the correct electron temperature. 
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Fig. 4. — An accretion-jet model for XTE J1550-564 at MJD 51,696. The dashed, thin solid 
and dot-dashed curves show the emission from the jet, the LHAF, and the truncated thin 
disk, respectively. The jet emission contributes significantly at the H, I, and V bands. On 
the other hand, the jet contribution is negligible in the X-ray band. This explains why we 
see secondary peaks at the H, I, and V bands with no corresponding peak at the X-ray 
band. See SJ3]for details. 



